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Abstract 

Inclusive jet production (e + e~ — > e + e~+jet+X) is studied in collisions of quasi-real 
photons radiated by the LEP beams at e + e~ centre-of-mass energies t/s^T from 189 to 
209 GeV. Jets are reconstructed using the k± jet algorithm. The inclusive differential 
cross-section is measured as a function of the jet transverse momentum, p!^ , in the range 
5 < p!£ < 40 GeV for pseudo-rapidities, rf et , in the range -1.5 < rf ct < 1.5. The results 
are compared to predictions of perturbative QCD in next-to-leading order in the strong 
coupling constant. 
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1 Introduction 



We have studied the inclusive production of jets in collisions of two quasi-real photons at 
e + e~ centre-of-mass energies t/s^ from 189 to 209 GeV, with a total integrated luminosity 
of 593 pb _1 collected by the OPAL detector at LEP. The transverse momentuuQ of the 
jets provides a scale that allows such processes to be calculated in perturbative QCD. 
Calculations at next-to-leading order (NLO) in the strong coupling constant, a s , for this 
process are available [U [2]. Comparisons of these calculations to the data presented in 
this paper provide tests of perturbative QCD for the observables and kinematical region 
considered. Leading order (LO) Monte Carlo (MC) generators are used to estimate the 
influence of soft underlying processes not included in the NLO calculation. Furthermore 
the measured jet cross-sections may be used in studies evaluating the hadronic structure 
of the photon, which are beyond the scope of this paper. 

The jets are reconstructed using the k± jet algorithm |3j. Inclusive jet cross-sections in 
two-photon interactions have previously been measured at ^fs^ from 130 to 136 GeV by 
OPAL [4J|, and at t/s^ from 189 to 209 GeV by L3 [5]. In [5] an excess at high transverse 
momenta was observed in the data over the QCD calculations, for kinematical conditions 
very similar to those used in the present paper. 

At LEP the photons are radiated by the beam electrons^], and carry mostly small 
negative four- momenta squared, Q 2 . In this paper events are considered only if the electrons 
are scattered at small angles and are not detected. Both photons are therefore quasi-real 
(Q 2 ~ GeV 2 ). The interactions can be modelled by assuming that each photon can 
either interact directly or appear resolved through its fluctuations into hadronic states. 
In leading order QCD this model leads to three different event classes for 77 interactions: 
direct, single-resolved and double- resolved, where resolved means that the incoming photon 
interacts through its partonic structure (quarks or gluons). 



2 The OPAL detector 

A detailed description of the OPAL detector can be found elsewhere [5]. The central 
tracking was performed inside a solenoidal magnet which provided a uniform axial mag- 
netic field of 0.435 T along the beam axis. Starting from the innermost components, 
the tracking system consisted of a high precision silicon micro-vertex detector, a preci- 
sion vertex drift chamber, a large volume jet chamber with 159 layers of axial anode 
wires and a set of z chambers measuring the track coordinates along the beam direction. 
The transverse momenta, p?, of tracks were measured with a precision parametrised by 

■"■OPAL uses a right-handed coordinate system where the z-axis points in the direction of the e _ beam 
and the x-axis points to the centre of the LEP ring. The polar angle 9 and the azimuthal angle <p are defined 
relative to the +z-axis and +x-axis, respectively. In cylindrical polar coordinates, the radial coordinate is 
denoted r. The transverse momentum is defined as the component of the momentum perpendicular to the 
z-axis. The pseudo-rapidity rj is defined as r\ = — lntan(0/2). 

2 Positrons are also referred to as electrons. 
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ct Pt /Pt = a/0-02 2 + (0.0015 • pt) 2 (pt in GeV) in the central region |cos#| < 0.73. The 
position of the primary vertex is determined from the tracks. 

The magnet was surrounded in the barrel region (| cos6*| < 0.82) by a lead-glass elec- 
tromagnetic calorimeter (ECAL) and a hadronic sampling calorimeter (HCAL), which in 
turn were surrounded by muon chambers. Similar layers of detectors were installed in the 
endcaps (0.82 < | cos#| < 0.98). The small-angle region from 47 to 140 mrad around the 
beam pipe on both sides of the interaction point was covered by the forward calorimeters 
(FD) and the region from 33 to 59 mrad by the silicon-tungsten luminometers (SW). The 
latter were used to determine the luminosity by counting small-angle Bhabha scattering 
events. 



3 Kinematics and MC simulation 

The properties of the two interacting photons (i = 1,2) are described by their negative 
squared four-momenta Q 2 and the invariant mass of the photon-photon system, W. Each 
Q 2 is related to the electron scattering angle 6[ relative to the beam axis by 

Q? = -fe-^) 2 ^2^(l-cos^), (1) 

where Pi and p\ are the four-momenta of the beam electrons and the scattered electrons, 
respectively, and and ~E\ are their energies. Events with one or both scattered electrons 
detected (single-tagged or double-tagged events) are excluded from the analysis. Driven 
by the angular acceptance of the FD and SW calorimeters, a value of Q 2 = 4.5 GeV 2 is 
used in this analysis as the maximum possible Q 2 . The median Q 2 resulting from this 
limit cannot be determined from data, since the scattered electrons are not tagged, but is 
predicted by MC simulations to be of the order of 10 -4 GeV 2 . The invariant mass of the 
photon-photon system, W, can be obtained from the energies and momenta (_Eh,Ph) of the 
final state hadrons. 

All signal and background Monte Carlo samples used for detector corrections and back- 
ground determinations were passed through a full simulation of the OPAL detector [7]. 
They are analysed using the same reconstruction algorithms as are applied to the data. 

The Monte Carlo generator PYTHIA 5.722 [H [9] was used to simulate the signal pro- 
cesses for the determination of detector corrections, as large samples with full detector 
simulation and reconstruction were available. For all other purposes the more modern 
PYTHIA 6.221 was used to generate signal samples. PYTHIA is based on LO QCD ma- 
trix elements for massless quarks with the addition of parton showers and hadronisation. 
The following generators were used for the simulation of the six background processes that 
contribute significantly after the event selection described below: PYTHIA for Z/7* — > qq 
and e+e" -> W + W"; BDK [10J for 77 -> r+r"; HERWIG p] for deep-inelastic electron- 
photon scattering (7*7); KORALZ [12J for Z -> t + t~ and GRC4F [TS] for e+e" -> e + e"qq 
events. 
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4 Jet definition and event selection 



The data presented were collected by the OPAL detector at centre-of-mass energies = 
189 - 209 GeV and represent a total integrated luminosity of 593 ptT 1 . For the purpose of 
this analysis, the difference between the data taken at the various values of ^fs^ is small 
and therefore the distributions for all energies have been added. The luminosity-weighted 
average centre-of-mass energy is 198.5 GeV. The efficiency to trigger jet events in the 
selected kinematical region is close to 100% [I]. 

In this analysis, a sum over all particles in the event or in a jet means a sum over 
two kinds of objects: tracks satisfying the quality cuts detailed below, and all calorimeter 
clusters, including those in the FD and SW calorimeters. A track is required to have 
a minimum transverse momentum of 120 MeV and at least 20 hits in the central jet 
chamber. The point of closest approach to the origin must have a distance of less than 
25 cm in z and a radial distance, do, of less than 2 cm from the z-axis. For tracks with a 
transverse momentum larger than 5 GeV, do is required to be less than 0.15 cm, to ensure 
a good momentum measurement. Calorimeter clusters have to pass an energy threshold of 
100 MeV in the barrel section or 250 MeV in the endcap section of the ECAL, 600 MeV 
for the barrel and endcap sections of the HCAL, 1 GeV for the FD, and 2 GeV for the 
SW. An algorithm is applied to avoid double-counting of particle momenta in the central 
tracking system and their energy deposits in the calorimeters [I]. The measured hadronic 
final state for each event consists of all objects thus defined. 

Events with at least one jet are first preselected before the final event selection based 
on maximum likelihood distribution functions [14] is applied. The preselection criteria are 
as follows. 



• Using the k± jet algorithm, the event must contain at least one jet with \rf et \ < 1.5 
and a transverse momentum pip > 5 GeV. In this algorithm the distance measure 
between any pair of objects {i, j} to be clustered is taken to be 

^ = min(^,p^)(4/^o) with R% = (A Vij ) 2 + (A<^-) 2 . (2) 

Throughout this analysis we set Rq = 1. The pr of the reconstructed jet, p!^ , is 
calculated as the sum of the transverse momenta of all the particles in the jet. 

• The total summed energy deposited in the ECAL and the HCAL has to be less 
than 80 GeV. This removes most of the hadronic Z decays, including events with a 
radiative return to the Z peak. 

• To remove events with scattered electrons in the FD or in the SW calorimeters, the 
total energy sum measured in the FD has to be less than 55 GeV, and the total 
energy sum measured in the SW calorimeter has to be less than 40 GeV. 

• The z position of the primary vertex is required to satisfy \z\ < 5 cm and the net 
charge Q of the event calculated from adding the charges of all tracks is required to 
be \Q\ < 5 to reduce background due to beam-gas interactions. 
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• To remove events originating from interactions between beam electrons and the beam- 
pipe the radial distance of the primary vertex from the beam axis has to be less than 
3 cm. 

• The event must lie in the allowed ranges for the input variables of the maximum 
likelihood selection, as detailed below. 

The final event selection uses seven input variables for the likelihood function: 

1. The visible invariant mass measured in the ECAL only, Wecax ( m t ne range 
[0-80] GeV); 

2. The visible invariant mass calculated from the entire hadronic final state, W Iec (in 
the range [0-120] GeV); 

3. The number of tracks (in the range [6-70]); 

4. The sum of all energy deposits in the ECAL (in the range [0-80] GeV); 

5. The sum of all energy deposits in the HCAL (in the range [0.1-55] GeV); 

6. The missing transverse momentum of the event calculated from the measured 
hadronic final state (in the range from zero to A /s^/2); 

7. To improve the rejection of background coming from hadronic Z decays, an invariant 
mass, Mjih2, is calculated from the jet with highest p 3 ^ in the event and the four- 
vector sum of all hadronic final state objects in the hemisphere opposite to the 
direction defined by this jet (considered in the range [0.1-100] GeV). 

In comparing the preselected events to MC simulations, the signal MC generator 
PYTHIA 5.722 underestimates the normalisation of the cross-section by about 50% in 
this process, and is scaled up accordingly. A similar deficiency was also observed in our 
previous study on di-jets [15]. Furthermore previous studies have shown that the prediction 
of MC generators for jet events in photon-photon collisions where one of the photons is 
virtual is too low by about a factor of two [16]. The prediction of the contribution from 7*7 
events has been scaled up accordingly, resulting in an adequate description of all quantities 
used in the analysis. 

Figured] (a) and (b) show two examples of the input distributions used for the likelihood 
selection, which is performed separately for events with the highest p 3 ^ smaller than 30 GeV 
and events with the highest p 3 ^ larger than 30 GeV. The region of high p 3 ^, , where most of 
the discrepancy with NLO QCD is observed in [5], is strongly affected by background from 
Z/7* — > qq which is not important at lower p 3 ^ . A separate optimization of the selection is 
hence necessary to maximize the reach of the analysis in p 3 ^ . The output of the likelihood 
functions for the data and all simulated processes is shown in Figured] (c) and (d). The 
cuts on the likelihood outputs are chosen to be 0.26 and 0.98 for the low and high p 3 ^ 
region, respectively. Applying these cuts reduces the background by 99.5% while reducing 
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the signal by 71% in the high p^* region; in the low p!^ region, these reductions are 91% 
and 27%, respectively. 

The p 1 ^ distribution after the event selection is shown in Figure EJ The dominant 
background at low p^ is due to 7*7 events, while for high p^ the background is dominated 
by Z/7* — > qq events. To measure the cross-section, the background is subtracted bin-by- 
bin. 

The measured transverse momentum distributions still have to be corrected for the 
losses due to event and track selection cuts, the acceptance and the resolution of the 
detector. This is done using the PYTHIA 5.722 signal Monte Carlo events which were 
processed by the full detector simulation and reconstruction chain. The data are corrected 
by multiplying the experimental distribution with correction factors which are calculated 
as the bin-by-bin ratio of the generated and the reconstructed Monte Carlo distributions. 
This method only yields reliable results if migration effects between bins due to the finite 
resolution of the measurement are small. The bins of the -distribution have therefore 
been chosen to be significantly larger than the detector resolution, obtained from the Monte 
Carlo simulation. 



5 Systematic uncertainties 

The systematic uncertainties for each bin in p^ can be found in Table [U The values 
for each bin were averaged with the results from its two neighbours (single neighbour 
for endpoints) to reduce the effect of bin-to-bin fluctuations. The sources of systematic 
uncertainties considered are given below. The total systematic uncertainty for each bin is 
calculated by adding the contributions in quadrature. 

• The absolute energy scale of the ECAL is known to about 3% [17] for the jet energy 
range in this analysis. To estimate the influence of this uncertainty, the energy scale 
of the data is varied by this amount, and the analysis is repeated. 

• To assess the uncertainty associated with the subtraction of background events, all 
backgrounds - except for 7*7 - have been varied by 10%. The prediction of the con- 
tribution from 7*7 events has been scaled up by a factor of two, as described earlier. 
By comparing the MC predictions in regions where this background dominates we 
conclude that this scaling factor can not be varied by more than about 30% while 
keeping a good description of the data. The scaling factor is varied accordingly. The 
resulting uncertainty is dominated by the distributions from 7*7 and Z/7* — > qq 
background. 

• To test the event selection's dependence on the simulation of the signal, the signal 
MC has been re-weighted to have a pO^ -slope in which it significantly either over- 
or underestimates the data at high pl^ and the analysis has been repeated. The 
difference between using the original MC and the re-weighted MC is included in the 
systematic uncertainty. 
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• The cut on the likelihood output value is varied down to 0.23 and up to 0.29 for the 
low p¥^ region and down to 0.88 for the high p!^ region. 

• The uncertainty on the determination of the integrated luminosity is less than 1%, 
and is neglected. 

6 Multiple parton interactions and hadronisation cor- 
rections 

The NLO calculations do not take into account the possibility of an underlying event, 
which leads to an increased energy flow and therefore to a larger cross-section above any 
given threshold in the jet transverse momentum. PYTHIA 6.221 has been used to study 
the effect of either considering (default) or leaving out multiple interactions for the signal 
MC. In PYTHIA the underlying event is modelled by multiple parton interactions (MIA). 
At the lowest transverse momenta considered the signal MC cross-section increases by up 
to 20% when including MIA. This effect reduces to less than 10% for transverse momenta 
larger than 7 GeV. 

The measured inclusive jet cross-section will be compared to NLO QCD calculations 
which describe jet cross-sections for partons, while the experimental cross-section is pre- 
sented for hadrons. There is as yet no rigorous way to use the MC generators to correct 
the NLO predictions for this process so that they can be compared to the data, because 
the partons in the MC generators and the partons in the NLO calculations are defined 
in different ways. But because the use of MC generators is the only available option so 
far, they are used to approximate the size of this hadronisation correction. Hadronisation 
corrections have been estimated with PYTHIA 6.22l£§ At p^ = 5 GeV the correction is 
about 15%. The correction decreases with increasing p^ and is below 5% in our study for 
p s ^ > 20 GeV. Disabling MIA in PYTHIA while determining the hadronisation corrections 
leads to values of the correction factors within 2% of those determined with MIA enabled. 



7 Differential cross-section 

Inclusive jet cross-sections have been measured for the photon-photon kinematical re- 
gion of invariant masses of the hadronic final state W > 5 GeV, and a photon virtual- 
ly Q 2 < 4.5 GeV 2 . The data are compared to predictions of PYTHIA 6.221 and NLO 
perturbative QCD [H[2]. 

The NLO cross-sections are calculated using the QCD partonic cross-sections in NLO 
for direct, single- and double-resolved processes, convoluted with the Weizsacker- Williams 

3 The use of the cluster fragmentation model as implemented in HERWIG as an alternative in studies 
carried out in [15) for similar kinematic conditions yields corrections compatible with or smaller than those 
determined using PYTHIA. 
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effective photon distribution. The hadronisation corrections discussed in the previous sec- 
tion are applied to the NLO calculation before it is compared to the data. The GRV-G HO 
parametrisation of the parton densities of the photon [18] is used with A^j S = 131 MeV. 
The renormalisation and factorisation scales in the calculation are set equal to . The 
cross-section calculations were repeated for the kinematic conditions of the present analy- 
sis. The calculations shown below are obtained from pQ. We have verified that using the 
independent calculation presented in [2] yields results within 5%, except in the lowest bin 
in pp*, where it predicts a cross-section about 25% higher. 

Figure [3] and Table [2] show the cross-section as a function of p!^ for \r]^ et \ < 1.5. Both 
PYTHIA 6.221 and the NLO calculation achieve a good description of the data, with the 
exception of the lowest bin in pl^, where the NLO calculation is too low. 

To facilitate a comparison with a recent measurement by the L3 collaboration, a mea- 
surement of the same quantity as presented in Figure [3] is shown in Figure H] for \rf et \ < 1.0. 
While the L3 data points are compatible with the present measurement, they lie below 
the OPAL data points at low p^ and above the OPAL data points at high pij^, and there 
is a discrepancy in shape between the L3 data and the NLO calculation. This difference 
in shape has been reported in the L3 publication and leads to a significant disagreement 
between the L3 data and the NLO calculation at the highest pip of up to 50 GeV stud- 
ied in [5]. The present analysis finds the region of > 40 GeV to be dominated by 
background and hence no measurement is presented for this region. 

In contrast to the conclusion in [5], the present analysis finds good agreement between 
data and calculations for p!^ of up to 40 GeV, leading to the conclusion that perturbative 
QCD in NLO is adequate to describe the process under study. 
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Ta&Ze 1: Systematic uncertainties on the inclusive jet cross-section in the individual p 3 ? 
bins for \rf et \ < 1.5. Values for \rf et \ < 1.0 are similar. 
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Table 2: Background fraction and inclusive jet cross-section for \r]^ et \ < 1.0 and \rj^ et \ < 1.5 
as a function of y 3 ^ ■ For the cross-section values the first uncertainty is statistical, the 
second is systematic. The uncertainty given for the background fraction is statistical only. 
The average value ofp 3 ^, (p 3 ?), is also given. 
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Figure 1: Example inputs to the likelihood functions: (a) shows the number of tracks in 
the event, and (b) the invariant mass of the system formed by the jet with the highest p 3 ^ in 
the event and the four-vector calculated from all objects in the opposite hemisphere as seen 
from this jet. Outputs of the likelihood functions: Plots (c) and (d) show the output of the 
likelihood functions for events with p 3 ^ < 30 GeV and p 3 ^ > 30 GeV, respectively. Events 
are selected with likelihood values larger than the cuts indicated by the arrows. The signal 
MC and the contribution of the 7*7 MC have been scaled up as described in section^ 
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Figure 2: Number of jets in eachp 3 ^ bin after the event selection. Monte Carlo distributions 
are normalised to the luminosity of the data. The statistical uncertainty of the data points 
is shown where larger than the marker size. 
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Figure 3: Inclusive jet differential cross-section, da/dp 3 ?, for all jets with \rf et \ < 1.5 
compared to NLO and PYTHIA 6.221 predictions. The hadronisation corrections described 
in section^ have been applied to the NLO results. The total of statistical and systematic 
uncertainties added in quadrature is shown where larger than the marker size. The inner 
error bars show the statistical uncertainties. The band on the NLO shows the uncertainty 
associated to the variation of the renormalisation and factorisation scale. 
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Figure 4 : Inclusive jet differential cross-section, da/dp 3 ^, for all jets with \rf et \ < 1.0 
compared to the results of the L3 collaboration, NLO and PYTHIA 6.221 predictions. The 
hadronisation corrections described in section^ have been applied to the NLO results. The 
total of statistical and systematic uncertainties added in quadrature is shown where larger 
than the marker size. The inner error bars show the statistical uncertainties. The band 
on the NLO shows the uncertainty associated to the variation of the renormalisation and 
factorisation scale. 
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